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A new approach to the automatic equalization of differentially phase- 
modulated data signals passing through a linear distorting medium is 
presented. The equalizer, which is of the transversal filter type, operates 
in the frequency passband and contains two sets of taps-in-phase and 
quadrature. A tap-rotation property of the equalizer is used to establish 
an absolute phase reference at the equalizer output, and once this crucial 
step is taken, the passbaird equalizer output can be used to automatically 
(as well as adaptively) adjust the tap weights so as to minimize a mean- 
square distortion function. The resulting algorithm requires correlating an 
error signal with the tap voltages, and thus it is possible to use a structure 
similar to that employed when equalizing baseband PAM. 

The generality of the approach makes it applicable to the equalization 
of any double-sideband data signal. 

I. INTRODUCTION 

Recently much attention has been given to high-speed synchronous 
data transmission via differential phase modulation and comparison 
detection. High-speed (above 2400 bits per second) transmission 
usually requires equalization 1 to compensate for the linear distortion 
introduced by the channel. This study will be concerned with the auto- 
matic and adaptive equalization of such data signals. At first glance, 
the nonlinear nature of the modulated signal would seem to preclude 
linear compensation. Upon closer examination it becomes apparent 
that a digitally phase-modulated signal, when resolved into in-phase 
and quadrature components, is linear in each component. Since this 
property is preserved after transmission through a linear medium, 
simultaneous linear equalization of each signal component becomes 
feasible. Due, however, to the purely quadratic nature of the channel- 
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dependent terms present in the differential detector output, a linear 
equalizer should 1 precede the detector-hence the equalizer operates in 
the frequency passband. 

Resolving the channel output into in-phase and quadrature com- 
ponents suggests that satisfactory performance will be obtained if the 
equalized in-phase signal is approximately Nyquist and the equalized 
quadrature signal has most of its samples close to zero. This is accom- 
plished by choosing an appropriate cost function (of the equalized 
samples) and using an appropriately structured equalizer capable of 
minimizing this function. We will focus our attention on a mean-square 
cost function and on a synchronous tapped delay line (TDL) equalizer. 
The equalizer has two branches-an in-phase and quadrature branch- 
which together perform the passband compensation. The delays are 
separated by a symbol interval, and each delay output is first multi- 
plied by a tap weight (each signal in the quadrature branch is also 
shifted by 90 degrees) and then added to the other delayed tap signals. 
The generality of the approach makes it applicable to the equalization 
of any double-sideband modulated data signal (e.g., combined ampli- 
tude and phase modulation, quadrature amplitude modulation). 

Though the differential detector output is independent of any 
absolute phase reference, it is convenient to select a cost function which 
depends on such a quantity. The awkwardness of this situation is re- 
solved by noting the tap-rotation property of the equalizer. Simply 
put, this means that any assumed reference phase angle can be 
"matched" by a rotation of the equalizer taps. This, in effect, permits 
the establishment of an arbitrary phase reference at the equalizer 
output. Equalization is accomplished with respect to this arbitrary 
reference while detection is done in an incoherent fashion. The tap 
adjustment algorithm is similar to that done in baseband PAM, i.e., 
tap signals are correlated with error signals. In fact, once the notion 
of a tap rotation is introduced, there is a convenient analogy to base- 
band PAM equalization. 

The basic system equations are indicated in Section II, and the 
equalizer cost function and structure arc described in Section III. 
Section IV describes methods for adjusting the equalizer taps when 
an ideal data reference sequence is available and also when the adjust- 
ments are made using random data. In Section V we consider the 
effects of frequency offset and phase jitter on the equalizer. 



t In order that the mean-square distortion be a quadratic, and thus easily mini- 
mized, function of the in-phase and quadrature pulse samples. 
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II. BASIC SYSTEM EQUATIONS 

In this section wo indicate the response of a linear passband channel 
to a phase-modulated data signal. We also make some observations 
concerning the choice of phase reference that will be useful in our 
discussion of passband equalization. 

The channel input is a phase-modulated (PM) data signal. This 
signal is given by 

Si(t) = Zp(t- nT) cos (bid + On) (la) 

= u in (t) cos bid — v in (t) sin m c t (lb) 

= r in (t) cos {bid + *i»(0)i (lc) 

where 6„ is the information symbol, 1/T is the symbol rate, u r is the 
carrier frequency, p(t) is the impulse response of the transmitter 
shaping filter,* and 

M in (0 ■ E p(t - nT) cos d n (Id) 

vUt) ■ E p(t - nT) sin 6„ (le) 

are respectively the in-phase and quadrature signal components. The 
envelope and phase of s,-(<) are given by 

r, a (0 = V (0 + y. 2 (0 (If) 

in in 

Vin(t) 

* in (t) = tan" 1 ——- (lg) 

Win(0 

By letting 

a n = cos0„ (lh) 

b n = sin0„, (li) 

we can write 

8<(i) = [E a n p(t - nT)] cos bid - [E b n p(t - nT)] sin w e t. (lj) 

n n 

From eq. (lj), we can see that the bandlimited signal s,-(i) is linear 
in both a„ and b,,* As we have previously remarked, it is this linear 



* We assume the usual Nyquist shaping in the sense that p(kT) = 0, k ^ 0. 

* Note that any double-sideband signal, such as quadrature amplitude modulation 
(QAM), can be written in the form of (lj). What we have to say in the sequel is 
sufficiently general to apply to any such double-sideband signals. 
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representation that makes linear compensation first plausible then 
possible. 

We note that the transmitter changes the value of 0„ once every T 
seconds, and that for differentially coded data, the customer infor- 
mation is the quantity 6„ — 0„-i. Applying s,(0 to a bandpass channel 
with impulse response 

2Fi(t) cos u> e t - 2F 2 (0 sin uj, 

produces an output signal 

s (t) = u (t) coswci — v (t) sin u c t, (2a) 

where 

u„(t) = fi(«)®«in(0 - F 2 (l)®v ia (t) 

v„(t) = F 2 (t)®u in (t) + F 1 (t)®Vin(t), 



(2b) 



with ©denoting the convolution operation. To obtain a more compact 
representation, we define the in-phase channel pulse 

x(t) =p(t)®F l (t), (3a) 

and the quadrature channel pulse 

y(t) = p(t)®F 2 (t). (3b) 

In terms of the data components 

a n = cos 8 n , 
b„ = sin 6„, 

and the above notation, the channel output is given by 

s o (0 = cos co c *[E a n x(t - nT) - b„y(t - nT)~\ 

n 

- sin »c«[L a»y(t ~ nT) + b„x(t - nT)~\. (4) 

The above in-phase and quadrature representation clearly indicates 
the effect of the channel on the transmitted signal. A "good" channel 
would be one which has a small quadrature pulse f (ideally all its 
samples would be zero) and an in-phase pulse which is Nyquist. In 
any event, it is important to note that in passband the channel dis- 
tortion appears in a linear fashion. 



t If the channel has even amplitude symmetry and odd phase symmetry about the 
carrier, then y(t) will be identically zero for all t. 
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Since the carrier phase reference is arbitrary, the output signal can 
also be written as 

s (t) = cos M + <*>)[£ a n x l *>(t - nT) - b n y^(t - nT)] 

n 

- sin (o> c f + 4>)[L a n y l *Kt - nT) + b n x^(t - nT)~\ (5) 

where ^ ( * } (0 and y {4,) {t) denote the in-phase and quadrature pulses 
with respect to the reference phase* <l>. Comparing the right-hand sides 
(RHS) of (4) and (5) gives 

r* ( * } (0"| r cos * sin **ir*(0"i 

Ly ( * ) (0J ~ L-sin* cos0_U(O-T 

Equation (6) indicates that changing the phase reference by radians 
corresponds to rotating the in-phase and quadrature signal vector 
[_x(t), y(t)2 by this amount. 1 We wish to emphasize the above in- 
terpretation, since we will want to consider the effect of such a rotation 
on several quantities of interest. We indicate the incoherent nature of 
the system by including an arbitrary phase angle when describing the 
channel output, i.e., the absolute phase of the received signal is un- 
known to the receiver. For incoherent reception, it is common to use 
differential coding of the data and differentially coherent (comparison) 
detection. We will assume this mode of transmission and detection. 

III. MEAN-SQUARE PASSBAND EQUALIZATION 

We wish to develop an adaptive equalizer which is capable of com- 
pensating for the linear (in-phase and quadrature) channel distortion, 
while keeping the equalizer adaption (settling) time small. "Mean- 
square" adaptive equalization 12 has been successfully used in systems 
using linear modulation to meet the above objectives. Our approach 
will be to extend, with appropriate modification, this technique to the 
problem at hand. After specifying the equalizer structure, we select 
an appropriate cost function and indicate how the equalizer processes 
the received random signal to adapt to the optimum configuration. 
We begin by considering the output of the differential detector. 



+ We reserve the right to suppress the phase reference when there is no possibility 
of confusion. Of course the phase reference can be chosen as convenience dictates. 
i Clearly the same relation also holds for the sampled values of x (t) and y (t) . 
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3.1 The Detector Output 

The equalizer is to be compatible with the standard comparison 
detector 1 which is shown in Fig. 1. It is straightforward to show that 
the detector outputs at the nth sampling instant are f 

I{nT +J,) = H COS (0,- - d,)[Xn-4Xn-j-\ + Vn-iVn-i-i} 
' i 

+ E E sin (0, - dj)[x n -,yn-j-i - y*-iX n -j-i] (7) 

» i 

Q(nT + t ) = E E COS (0, - 9j)\jfn-iXn-i-l ~ Xn-»1/n-i-l] 

I J 

+ E E sin (0,- - 0j)[Xn-iXn-i-i + y«_,y«_/_i] (8) 

i i 

where 

x(nT + t„) = x n 
and 

y(nT + t„) = y„. 

We note that the intersymbol interference at the detector output is of 
a quadratic nature and thus is quite different than that encountered 
in linear modulation. It should also be noted that the channel-depen- 
dent terms in (7) are inner products, and hence are independent of the 
phase reference. The terms in (8) are components of cross products 
and hence also independent of phase. This is a manifestation of the 
differential nature of the receiver, i.e., only quantities which are in- 
variant with respect to a change in phase reference can affect the 
detector output. 

If we are fortunate enough to have a distortionless channel, i.e., for 
all integers n, 

x„ = 5 no (9a) 

y n = 0, (9b) 

then detection can be simply accomplished by noting that 

I(nT + to) = a n a n -! + &„&„_! = cos(0„ - 0„_i) (10a) 

Q(nT + h) = b n a„-i - a„6„_i = sin(0 n - 0„-i) . (10b) 

Unfortunately, (9) will not be satisfied for an arbitrary channel. As 
mentioned above, the function of the proposed equalizer is to linearly 



t We have, of course, neglected noise and have suppressed the phase reference. 
We also wish to emphasize that a desirable sampling epoch t a should be determined. 
As usual, this is a problem in its own right, but is not considered further in this paper. 
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Fig. 1 — Comparison (differential) detector. 

operate on the in-phase and quadrature samples (the x„'s and y n 's) 
to produce an overall system response the samples of which (denoted 
by g n 's and h„'s) approximately satisfy eq. (9). The quality of this 
approximation is measured by the value of an appropriate cost func- 
tion. We will have more to say about this later. 



3.2 Equalizer Structure 

Since we would like to equalize both the in-phase and quadrature 
channel samples, a structure which combines two signals in quadrature 
is suggested. The equalizer, which is shown in Fig. 2, consists of a 
tapped delay line (TDL) where each of the 2N + 1 delay outputs is 
fed into two branches. In each branch the delay outputs are multiplied 
by tap weights (c, and di) and then summed. The equalizer output 
q(t) is the sum of the upper output qi(t) and the Hilbert transform of 
the lower output q->(t). The 90-degree phase shift at the output of the 
lower branch provides the quadrature signal. Thus the taps are to be 
adjusted to minimize an appropriate (cost) function of the equalized 
in-phase and quadrature samples. We remark that a linear equalizer 
placed in the passband can compensate for linear channel distortion, 
while a linear equalizer placed after the comparison detector cannot 
hope to compensate for the quadratic distortion present in the detector 
output. Thus the equalizer precedes the detector and operates in the 
frequency passband. We now describe the equalizer operation in detail. 
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Fig. 2 — The passband equalizer. 

We can write the equalizer output, q(t), as f 
q(t) = coso> c *[£ a„g(t - nT) - £ b n h(t - nT)~\ 

n n 

- sin uJCZ a„h(t -nT) + £ Kg(t - nT)~\ (11) 



where the in-phase and quadrature (equalized) signals, g(t) and h(t), 
are given by 1 

g(t) - £ c n x(t -nT) - £ d n? /« - nT) (12a) 

n=-AT H H 

h(t) = £ d n x(t -nT) + £ CfK* - nT). (12b) 

n=-AT n=-tf 

We emphasize that s,(0 and q(t) are denned with respect to the same 

+ Assuming that u c and T are chosen such that to e T is an integer multiple of 2ir. 
If for some reason this is not convenient (perhaps a particular carrier frequency is 
desired), then u c T can be chosen to be any convenient angle. Some additional book- 
keeping is then required at the receiver. 

* The time reference has been taken at the center tap of the equalizer. 
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phase reference, but that the phase reference has not been indicated 
in the above equations. For a sampling epoch t , we denote the equal- 
ized samples by 



Qi-gUT + Q, hj = hUT + t ). 



(13) 



Thus with the reference phase taken to be <f>, the equalized samples 
are related to the channel samples and tap weights by 



.v 



N 



(#) ^ (*) ^ . (*) 

(Jj = L CnXj-„ - 2- dnVj-n 

-N -N 

<*) * (*) * <*) 

hj = 2- dnXj-n + 2w CnVj-n. 

-N -N 



(14a) 



(14b) 



Again we reserve the right to suppress the representation phase in 
writing (14). We now make some observations which indicate the 
effects of a change in reference phase and of a "tap-rotation" on the 
equalized samples. 

Denoting the channel phase reference by <f>, we can write the equalizer 
output as 

q{t) = cos {u e t + 0)[D a„gW(t - nT) - £ b n h^\t - nT)~] 

- sin (« e * + 0)[L a„hi*>(t - nT) + £ &„</<♦>(< - nT)~], (15) 



and it is easy to see that 
p/<*>(0l 



cos<£ sin <P~\\~g(t)' 
_ — sin <j> cos0_ _h(t)_ 



(16) 



Thus we again observe that changing the reference phase by <f> radians 
corresponds to rotating the in-phase and quadrature signals (and 
samples) by this amount, i.e., each two-tuple (g ( f, h l f) is rotated by 
radians. Suppose that while keeping the phase reference fixed, we 
transform the taps (c, and (/,) to new values (c, and d f ) via the rotation 



- [ 



cos 



sin 



— sin 9 cos 0_ _di_ 



i= ~N, ...,0, ---N. (17) 



We can now express the equalized samples g, and hi (corresponding 
to Zi and di), using (14) and (17), as 



[a 



cos 6 sin 
— sin 6 cos 6 



ia 



i = - at, ...,0, ---N. (18) 
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Thus rotating the equalizer taps by will rotate the equalized samples 
by the same amount. f It is important to note that the original and 
(tap) rotated equalized samples are defined with respect to the same 
phase reference (0). This is in contrast to the rotation introduced 
when the reference phase is changed. We will use these notions in 
describing the equalizer adjustment algorithm, in fact, the key to the 
simplicity of the equalization system lies in recognizing that the 
possibility of tap-rotation allows us to fix the reference phase and 
still keep a rotational degree of freedom between all (</,, hi) pairs. The 
reference phase will be suppressed in the sequel. 

It remains to choose an appropriate cost function (of the equalized 
samples) and to describe an iterative procedure for determining the 
tap weights that minimize this function. 

3.3 The Cost Function 

The cost function should be such that after minimization, the 
equalized channel is approximately ideal (i.e., the intersymbol inter- 
ference has in some sense been suppressed). There are, of course, 
many such functions-thus our selection is influenced by additional 
considerations. 

We begin by noting that one possible task to set for the equalizer 
would be the minimization of 

D = Z gl + L hi (19) 

subject to appropriate constraint on g and ho. 1 We observe that since 
D is independent of an equalizer tap rotation, the minimum value of 
D is the same for any constraint that satisfies 

g\ + hi = 1, (20) 

including the choices 

go =1, h = 0. (21) 

Clearly a choice such as (21) fixes the reference phase used to describe 
the equalizer output, and a "rotation" of the equalizer taps can assure 
that it is met. Ultimately, the information needed to adjust the 
equalizer taps will indeed be obtained with respect to a particular 

t It should be clear that if there were no d tap, we could not perform this rotation. 

* We note that this cost function would also be appropriate for a data signal which 
employs combined AM and PM. In fact, it is a simple matter to extend all of our 
results (i.e., the equalizer structure and algorithm) to this signaling format. 
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phase reference [in fact, we use the one implied by (21)]. Hence the 
tap adjustments will be made using coherently obtained information. 

An additional interpretation of (21) is that the tap vectors are 
always required to lie on the hyperplanes g = 1, h = 0. Unfortu- 
nately, minimizing a cost function subject to this (hard) constraint 
cannot be done in real-time. This is due to the awkward requirement 
that the tap vectors be adjusted while constrained to lie on hyper- 
planes. For example, suppose we were using a steepest descent 3 
algorithm to adjust the taps, then we would need the projected (onto 
the constraint hyperplanes) gradient of D at every iteration. 1 At 
present it is not known how to easily generate such a gradient from 
the available circuit voltages. 

We circumvent this difficulty by imposing a quadratic penalty (a 
soft constraint) when the variables do not satisfy the constraints. 3 
Thus the equalizer will try to minimize the unconstrained mean-square 
distortion 

8 - E g\ + E h\ + (1 - g ) + hi (22) 

It should be noted that for small distortion, the optimum D [subject 
to (21)] and 8 will be essentially identical. 1 

A useful consequence of the notions of phase-reference rotation and 
tap rotation will now be demonstrated. Suppose the phase reference 
has been chosen and is fixed at this value. Then 8, which can be re- 
written as 5 

6 = E (9* + *!) + 1 - 2<7o, (23) 

nil i 

clearly depends upon the chosen phase reference, since g does. How- 
ever, since we have shown that a simple tap rotation can be used to 
effect any desired phase, we are guaranteed that the minimization of 
8 will be over all reference phase angles. That is, while (23) is in general 
not reference-phase invariant, the minimum of (23) over all tap settings 
is reference-phase invariant. A further manifestation of the tap- 
rotation property is the following necessary condition for the optimum 
tap setting : The optimum tap setting is such that h = 0. To see that 
this must be the case, suppose that the taps have settled down and 

+ A more detailed discussion of this point is available in Ref. 4. 
* See Ref. 4 for a discussion of this type of subject for baseband PAM. 
§ Since the phase reference has been fixed, we do not indicate this quantity when 
writing the equalized samples. 
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Fig. 3 — Illustration of the tap-rotation property. 

that h t* 0. Since E.fl? + ft? is unchanged under a tap rotation, the 
equalizer taps could now be rotated so as to further minimize & by 
maximizing g (i.e., making ho = 0). The rotation of the vector (9-0, ho) 
into the vector (^gl + h 2 Q , 0) is shown in Fig. 3; thus the required 

rotation is tan -1 ( W<7o) • 

Recalling that a differential detector is only affected by phase- 
invariant quantities, we observe that if it were not possible to perform 
the above sort of minimization, then the utility of the structure and 
cost function would be suspect. 



IV. EQUALIZER TAP ADJUSTMENT 

Now that we have selected a cost function, we wish to demonstrate 
how the equalizer taps should be adjusted to minimize this function. 
We proceed in a manner reminiscent of mean-square equalization for 
baseband PAM-in fact it is shown that the tap increments are ob- 
tained by cross correlating error signals with delayed equalizer inputs. 

4.1 Deterministic Aspects 

We first wish to show that our criteria, S, is an easily minimized 
function of the tap weights. We denote the tap vectors by 



(24) 



c T = (c- N , • • •, Co, ••', Cn) 

d T = (d- N , ■ ■ ■, do, ■••, dy) 

where the superscript denotes transpose. It is convenient to introduce 
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the channel correlation matrices 1 

L-^Jij = 2^ Xn—iXn—j 
n 

[Y~\ii - £ Vn-iVn-j 

[■£]« = E (Vn-iXn-j - Xn-xUn-j), 
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(25a) 
(25b) 
(25c) 



and the truncated channel vectors 

x T - (.r A -, • • -Xo, ■ ■ -X-n) 

y T = (y-N, ■•■yo, ■■•vn). 

We can now write the mean-square distortion, in matrix notation as 



(26) 



8 = (c T , d T ) 



X + Y 



K 



_ - K X + Y_\d 



CX -2(x^, -y T )(°\ -i l. i-.!7) 



A few comments are in order on the above representation. We first 
note that X and Y are the in-phase and quadrature channel correlation 
matrices. Thus if the quadrature component is zero, we need only the 
c taps. The cross-correlation between the in-phase and quadrature 
channel is measured by the skew-symmetric matrix K. By introducing 
the bandpass channel correlation matrix 



A = 



■X + Y K 
L - K X + Y_\ ' 

the augmented channel vector 

Z T = (x T j _ yT)j 

and the composite tap vector 

b T = (cW, 
the mean-square distortion can be expressed as 
S = b T Ab - 2b T z + 1. 



(28) 

(29) 
(30) 
(31) 



Conceptually at least, the optimum tap settings are given by setting 
the gradient of 8, with respect to b, to zero. This gives 1 



Ah - % = 



(32a) 



! C-Y],j denotes the ijth element of the matrix X. 
x Posith 



tive definiteness of -4 is assumed. For the optimum tap setting it is easy to 
see that 8=1- z T A~ l z. 
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or 

b* = A~ x z (32b) 

as the desired tap settings. 

If a gradient (steepest-descent) algorithm is used to minimize 1 6 
by iteratively adjusting b, then the speed of convergence of the 
algorithm (also referred to as the equalizer settling time) is determined 
by the matrix A. A simple bound on the speed of convergence is given 
in Ref. 4. 

1 / AminN" „ 

SnS (l -" ) So, 

A m in \ Aniax/ 

where X max and X min denote, respectively, the maximum and minimum 
eigenvalues of A, and 8 n is the mean-square error after the ?ith iter- 
ation. The ratio p = (XmaxAmin) was computed for several voice-grade 
telephone channels (note that the closer p is to unity, the smaller the 
bound on settling time). Typically, the effect of the quadrature channel 
was small and the A matrix was diagonally dominant with p in the 
range from 2 to 8. Chang 5 has computed p for baseband PAM trans- 
mission, and comparing our numerical values with his indicates that 
for most channels the passband equalizer will settle rapidly. In other 
words, there are no special phenomena arising in the equalization of a 
bandpass channel which might lead one to expect a larger settling time 
than that observed for baseband channels. Though (32b) tells us what 
the optimum settings should be, we have not as yet indicated how 
these settings can be generated from available circuit voltages. 

4.2 Real-Time Tap Adjustments Using an Ideal Reference 

Ultimately, the cost function actually minimized is determined by 
the ease with which the optimum tap settings can be obtained in real 
time. We first note that if, as assumed, A is positive definite, then 8 
is a convex function of the tap vector b. Since, as is well known, a 
convex function has a single minimum, a gradient (steepest-descent) 
algorithm can be used to adjust the taps. If h w denotes the tap vector 
after the fcth adjustment, then the next tap setting is given by 

b<*+D = b (t) - a (fc) V£(b (t) ). (33) 

where a {k) is a gain (or step size) and V 8 is the gradient of 8. Note 



t The details of how one uses such an algorithm in this context are described in the 
next section. The reader not familiar with this type of algorithm should read Section 
4.2 before reading the present discussion. 
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that the algorithm "turns itself off" only when the gradient is zero. 
For a careful choice of the sequence {a lk) \, the above algorithm will 
converge to the optimum tap setting. Clearly the major problem 
associated with a real-time implementation of (33) is generating the 
gradient-or a good (statistical) approximation. To obtain an estimate 
of the gradient, we parallel the approach taken in baseband PAM. 
The gradient of & can be rewritten in component form, as 

dS 

— = E (fJnXn-i + knVn-i) ~ X-i (34a) 

dCi 

r)& 

— = E (h n xn-i - gnijn-i) + y-i i = - N, ■ • -0, • • -N. (34b) 

Bdi 

We remark that a reasonable statistical approximation (i.e., an esti- 
mate) of the true gradient would be one whose average value is the 
right-hand side of (34). To do this, we first use the sampled equal- 
izer output and an ideal data reference to generate a quantity, the 
average value of which is precisely our cost function. Once we have 
done this, we interchange the (linear) operations of ensemble averaging 
and differentiation to obtain the desired estimate. 
The equalizer output samples are given by + 

q(jT + t ) b q s = cos (u e t + 0)[E {a*Qi-* ~ &»fy-n)] 

- sin (u c to + 0)[E (a*hj- n + &n0;-n)], (35) 
and for perfect equalization these samples would be 

q, = cij cos (o) e to + ^) — bj sin (co c t o + <£) = cos (u) c to + <t>j + <^>) • (36a) 

If we have an ideal reference (i.e., a, and bj are known at the receiver), 
then since the tap-rotation property permits the equalizer to obtain 
any phase, we have the liberty to choose <t>. In other words, we estab- 
lish a phase reference at the equalizer output, and we force (via the 
cost function penalty) the equalizer to find the phase angle of the 
vector (g { o\ h$) which minimizes the cost function. The equalizer 
does this by rotating the taps until h l o = 0. For convenience we choose 
the reference phase, <f>, so that 

qj = (a y - bj), (36b) 



1 Again, since the reference phase is arbitrary, we are free to choose as convenience 
dictates. 
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i.e., w c to + has been taken to be 7r/4. Note that we do not need an 
oscillator to generate the ideal reference. By analogy with baseband 
PAAI, we next consider the mean-squared difference between the 
sampled equalizer output (g,-) and the ideal (or desired) sequence 
(g*), i.e., EUq, - g*) 2 ]. With 

a 2 = E[al~\ = E[bl], (37a) 

the independence of the data sequences is used to show that 

*L<Ml ~ <?*) 2 ] = «*CE fo» + *!) + 1 - 2 ? o] = a&. (37b) 

n 

This is an extremely useful observation, since we now have available 
a quantity (g,- — g*) 2 , the average value of which is proportional to the 
cost function. Interchanging the linear operations of expectation and 
differentiation allows us to write 



d .2 T * tyi 

— £[(<?, " Qj) 1 = 2E\ ( Qj - qj) — 
dCi L dC(_ 

and 

— E£(qj - qj) ] = 2E\ (qj - q,) — 
ddi L ddi_ 



(38a) 
(38b) 



Note that q } - — q* is the error (at the equalizer output) at the jth 
sampling instant. It is easy to see that the differentiated terms can be 
written as 

dqj 

= COS (o) c t + <£)[X a n Xj-i~n — L hnVj-i-n'] 

dCi n n 

— Sin (u c to + (/»)[E dnljj-i-n + L b n X,-i-nl i ( 39a ) 

n n 

and 
dqj 

— - = - COS (w c to + 4»)[L dnVj-i-n + L bnXj-i-nl 

ddi n 

— sin {o) c t Q + </>)[5I a,nXj-i-n — 2 b„yj-i- n ~]. (39b) 

n n 

We recognize that dq,/dci is the jth channel output sample delayed by 
%T seconds, and thus is available at the z'th delay element of the 
equalizer. We also note that dqj/ ddi is precisely dqj/dCi, but with the 
carrier phase delayed by -w/2 radians, and thus is available at the out- 
put of the ith. delay followed by a tt/2 phase shifter. The estimated 
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gradient, which is given by the bracketed terms on the RHS of (38), 
can thus be obtained by correlating (muliplying and perhaps then 
averaging) the error signal (qj — q*) and the tap outputs. 1 As the 
equalization becomes better, the ideal reference can be replaced by 
decisions, i.e., the receiver operates in a decision-directed mode. This 
then gives a complete description of the passband equalization of the 
PM data signal when an ideal reference is available. 

At this juncture we wish to make the following remark. Suppose the 
channel is initially perfect [in the sense of eq. (9)3, but the phase 
reference is such that h is not zero:* What does the equalizer do? We 
first observe that the equalizer will rotate the taps until h is zero, 
and thus lock in on the established reference phase. However, the 
differential detector output, which is insensitive to such a rotation, 
will have been error free from the outset. 

4.."J A Decision-Directed Tap Adjustment Procedure 

Suppose an ideal reference has been used to start-up the equalizer. 
It is desired that the equalizer then be capable of using decisions in 
place of an ideal reference. This mode of operation has been dubbed 
adaptive or decision-directed equalization. 1 The correction term needed 
for the tap adjustment algorithm is the product of the equalizer error 
signal (qj — q*) and a delayed equalizer input. For decision-directed 
operation, q* (which is a, — b,) would be replaced by the decided value 
of Qj. For four-phase operation, a, and bj are binary, thus a 7 — 6 ; 
= cos(0> + w/4) takes on the values —0.707 and 0.707. Thus g> can 
be obtained as the output of a threshold device (with a threshold at 
zero). 

V. FREQUENCY OFFSET AND PHASE JITTER 

In this section we consider the effect of two common transmission 
impairments, frequency offset and phase jitter, on the operation of the 
equalizer. 

5.1 Frequency Offset 

Frequently, the transmission media perturbs the received carrier 
frequency w t . so that it differs from the transmitted carrier frequency 

* Remembering that to obtain dqj/ddi the ith tap signal has been shifted by w/2 
radians. 

1 This could come about, for example, by a phase-hit impinging on an otherwise 
ideal channel. 
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2w/T by A Hz. f As a result of this frequency offset, each tap signal has 
a different carrier phase. To obtain a more detailed understanding of 
this effect on the equalizer, we introduce the following rotated equalized 

on TT\ Tl IpQ 

|~0»(A)1 [" cosnA sinnA-p„n 

l_MA)J L — sin n A cosnAjL/inJ 

and . 

pm.n(A)-| r cosmA sin mAnp n (A)-| 

U m ,„(A)J L-sinwA cosmAjUn(A)J ' 

where g n and h n are the equalized samples in the absence of frequency 
offset. 

It is straightforward to show that the equalizer output, at the mth 
sampling instant, can be written in terms of the above quantities as 

q{mT) = E Om-«I m .»(A) - £ & m -JL, n (A). (42) 

n n 

The equalized samples [which are precisely the | m .„(A)'s and the 
A m , n (A)'s] are thus a cascade of these two rotations. The rotation 
described by (40) is time-invariant and amounts to presenting the 
equalizer with a slightly different channel which is obtained by rotating 
the vector (g„, h n ) by nA degrees. By observing that 

gl(A) + hl(A) = g\ + hi 
and , , y* 6 ) 

O (A) = 0o, 

it is clear that this effect can be neglected since the cost function is 
invariant under such a transformation. 

The rotation described by (41) is time-varying and indicates that at 
the mth sampling instant, each channel pair (g n (A),h„(A)) is rotated 
by mA radians. Since each channel pair is rotated by the same amount, 
an equalizer tap-rotation can compensate for this rotation, the only 
requirement being that the equalizer settling time be much smaller 
than the period of rotation 1/A. Thus if the settling time of the equal- 
izer is small, the equalizer automatically tracks the frequency offset 
by rotating each tap pair (c,,d,) at the offset frequency A. 

5.2 Phase Jitter 

Phase jitter is an additive random component 6(t) often present in 
the phase angle of the channel output, and is characterized as a low- 

t Typically A is less than 0.5 Hz. 
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pass process with energy up to 150 Hz. In this section we give only a 
suggestive description of the effect of phase jitter on the equalizer 
operation, since a detailed study appears to be difficult and could well 
be the subject of a separate investigation. If one were to assume the 
jitter to be constant along the length of the equalizer, i.e., 

e(mT - if) = 9(mT) i= - N, • • •, 0, ■•-, N, 

then at the ?nth sampling instant, each equalized sample pair (g n , h n ) 
is rotated by 6{mT) radians. Such an assumption could only be valid 
for the extreme low-frequency components of the phase jitter, of 
course, and these would be successfully "tracked" by slow to-and-fro 
motions of the tap settings. The higher-frequency components could 
not be followed by the equalizer, but might be modeled by noise sources 
at each tap causing fluctuations of the individual settings about the 
optimum values. In such a model, it would be realistic to include cor- 
relations among these fluctuations at each spectral component and one 
would expect the correlations to decrease with increasing frequency. 

VI. SUMMARY AND CONCLUSIONS 

An approach to the mean-square equalization of a differentially 
phase-modulated data signal has been presented. An equalizer, of the 
transversal filter type, has been proposed which operates in the fre- 
quency passband and contains two sets of taps-in-phase and quad- 
rature branches. By exploiting the tap-rotation property of the equal- 
izer, a phase reference is established at the output of the equalizer. 
Among the manifestations of the tap-rotation property are the ability 
to control the equalizer by using coherently obtained output samples 
and the ability to track small amounts of frequency offset and phase 
jitter. The equalized output is used to automatically (as well as 
adaptively) adjust the equalizer taps so as to minimize a mean- 
square distortion function. The required operations (correlating an 
error signal with tap voltages) are those performed when equalizing 
baseband PAM. Thus much of the existing knowledge concerning the 
technology of baseband PAM equalization can be applied to the 
equalization of phase-modulated data signals. 

We make two final remarks. The equalizer structure and tap adjust- 
ment algorithm can be applied, with minor modification, to the 
equalization of any double-sideband modulated data signals, e.g., 
combined amplitude and phase modulation. As is the case with all 
equalization systems, the precise dynamic behavior of the proposed 
equalizer can only be studied, at present, by experiment. 
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